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A  Parameterized  Procedure 
for  Determining  Real  Height  From  lonograms 
by  Use  of  Generalized  Parabolic  Profiles 


1.  INTRODUCTION 

1  In  previous  reports  ’  an  analytic  procedure  was  developed  for  determining 
real  heights  from  ionograms  by  directly  solving  the  integral  relation  between  real 
height  and  virtual  height.  The  real  height  could  then  be  obtained  by  a  direct  inte¬ 
gration  of  the  virtual  height  over  the  range  of  frequencies  encountered.  In  a  later 

-3^ — 

report  the  method  was  applied  to  typical  ionograms;  the  required  integration 
being  performed  numerically.  In  addition,  a  simplified  procedure  was  developed 
in  which  an  ionogram  could  be  fitted  to  a  parametric  series  of  polynomials  and  the 
integration  required  for  dete rmination  of  real  height  performed  analytically. 

-  In  the  pr.oD^rrt  report  generalizes^this  procedure  by  modeling  electron  density 
profiles  in  terms  of  a  parameterized  series  of  generalized  parabolic  profiles. 

The  profiles  are  then  converted  to  virtual  height  or  ionogram  form.  The  real 
height  corresponding  to  a  given  ionogram  is  obtained  by  comparison  of  the  ionogram 
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2)  with  the  parameterized  profiles.  Since  the  real  height  for  a  given  parameterized 
profile  is  related  in  a  very  simple  way  to  its  virtual  height,  the  integration  process 
required  for  the  previous  polynominal  method  is  hereby  avoided.  In  addition, 
since  the  parabolic  profile  has  a  vertical  slope  at  the  maximum  electron  density, 
it  provides  a  better  fit  to  ionospheric  data  than  a  polynomial.  tJ 


DETERM1P 


l ABO LIC  PROFILES  J 


ETERMINATION  OF  VIRTUAL  HEIGHT  FOR  PARABOLIC 


The  parabolic  profile  frequently  utilized  for  ionospheric  analysis  is  of  the  form 


!- JL 

NP 


where  N  is  the  electron  density,  Np  the  maximum  value  of  N,  hm  the  maximum 
height  of  the  profile,  and  hQ  the  initial  height  at  which  N  goes  to  zero.  This  pro¬ 
file  is  characterized  by  a  small  but  non-zero  slope  at  hQ  and  a  vertical  slope  at 
Np.  It  is  convenient  to  write  Eq.  (1)  in  terms  of  frequencies  by 


where  fpj  is  the  plasma  frequency,  f  its  maximum  value,  X  is  the  ratio  f^j/f^,  b 
the  ratio  f/fp,  and  1  is  the  probing  frequency.  The  foregoing  profile  is  character¬ 
ized  by  a  very  small  rate  of  increase  of  height  (and  hence  of  virtual  height)  at  low 
frequencies.  We  shall  refer  to  this  variation  of  f^  as  a  f^  profile.  In  order  to 
cover  ionograms  characterized  by  a  more  rapid  increase  of  height  at  low  fre¬ 
quencies,  it  is  useful  to  also  consider  the  profile  characterized  by  a  linear  varia¬ 
tion  of  frequency  (referred  to  as  a  linear  fN  or  f^  profile). 
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We  now  generalize  the  profiles  given  by  Eqs.  (2)  and  (3)  by  replacing  the  power  2 
by  a  parameter  n.  As  n  increases  the  slope  of  the  profile  decreases.  We  have 
found  that  a  range  of  n  from  1  to  4  is  generally  sufficient  for  the  ionograms 
encountered  in  practice.  Accordingly,  Eqs.  (2)  and  (3)  are  replaced  by 


*/«*  A'  V*  A’ 


(4) 


(hm  -  h)  =  a(l  -  b  X) 


(hm  -  h)n  --  a(l  -  bX1/2>  . 


(5) 


where  a  -  (hm  -  hQ).  A  plot  of  Eq.  (4)  exhibiting  the  form  of  the  generalized 


parabolic  curves  for  several  values  of  n  is  given  in  Figure  1.  For  subsequent  use 
we  write  Eqs.  (4)  and  (5)  in  the  form 


h„  “  h  -  a[  1  -  (1  -  b2)1/"] 


l1/0! 


hv  "  h0  =  afJ  “  (1  -  b)  1  • 


(6) 
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Figure  1.  Plot  of  Generalized  Parabolic  Curves  for 
Several  Values  of  Parameter  n.  The  electron  density 
N  has  been  normalized  to  its  maximum  value  N 


The  virtual  height  h * ( f )  corresponding  to  a  density  distribution  X  is  given  by 


",<ri  -  if  m»  ■  at  <»» 


where 


J  (1  -  xW2 


h  is  the  reflection  height,  and  a  is  a  parameter  incorporating  the  effect  of  the 
magnetic  field.  A  plot  of  <>  as  a  function  of  frequency  for  several  values  of  the 
angle  o  between  the  ray  and  the  nagnetic  field  is  given  in  Reference  3,  Figure  3. 
For  the  profiles  defined  by  Kqs.  (4)  and  (5),  we  may  write  Eq.  (9)  in  the  form 
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o  (1  -  b  X) 


ab  r  ( 1 

n  J  (1  - 


2  °!2 
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where  u  f^/f,  and  we  now  distinguish  between  the  f^  and  f^  cases  by  use  of  I 
and  J.  In  gener  al,  the  integrals  in  Kqs.  (10)  and  (11)  are  not  analytically  tract¬ 
able  for  arbitrary  n,  but  may  be  approximated  by  integrals  that  can  be  evaluated. 
We  first  write  Kqs.  (10)  and  (11)  in  the  form 
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-  XP  n (1  -  X)  n 
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(1  -  bu) 


and  note  that  the  quantities  (1  -  X)^2^  ^+n,  (1  -  u)(°/2)-l+n.  an(j  (j  +  / 2 

are  slowly  varying  functions  over  the  r  ange  of  integration  and  may  be  replaced  by 
average  values.  We  accordingly  approximate  Kqs.  (12)  and  (13)  by 
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where 
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A  numerical  check  of  the  exact  and  approximate  integrals  for  several  values  cf  b 
shows  that  the  approximate  forms  are  within  a  few  percent  of  the  corresponding 
exact  integrals. 

Carrying  through  the  differentiation  indicated  in  Eq.  (8),  the  virtual  height 
h'(f)  takes  the  form, 
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Since  we  shall  be  taking  rational  values  of  n,  the  power  1  -  1/n  may  be 
written  as  p/q,  where  p  and  q  are  integers.  Rationalizing  the  integrands  in  I  and 
J  by  use  of 


transforms  the  I  and  J  integrals  to  the  forms, 
2  2  f  VP+C!“1 

fN  •  Ji  =  <i  -  b  >  q  f  ?  dy 
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which  can  easily  be  integrated. 


3.  EVALUATION  OF  VIRTUAL  HEIGHT  FOR  SQUARE  fN  PROFILE 


3.1  n  =  4 

The  I  integrals  and  h'(f)  are  given  by 


I  =  4(1  -  b2)  f  - y-TT  dy 

1  o  (1  ‘  b  Vl 


4  =  4  /  - TT  dy 
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h'(f)  =  |  ab2g(I1  +  12)  . 


which  upon  evaluating  1^  and  Ig  yields 


h'(f)  =  rag  + 
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For  small  b,  Eq.  (34)  has  the  series  form 
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For  this  case  we  obtain 
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For  this  case,  the  integral  I  in  Eq.  (10)  can  be  evaluated  directly  to  yield. 
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which  gives  a  virtual  height 
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4.  EVALUATION  OF  VIRTUAL  HEIGHT  FOR  LINEAR  fN  PROFILE 
4.1  n  =  4 

The  J  integrals  and  h'(f)  are  given  by  (note  that  the  J  integrals  are  replicas 
2 

of  the  I  integrals  with  b  replaced  by  b) 
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The  integral  J  in  Eq.  (9)  can  here  be  evaluated  exactly  to  yield 
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and  a  virtual  height 
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5.  APPROXIMATE  SOLUTION  FOR  n  CLOSE  TO  UNITY 


The  integrations  required  in  the  foregoing  analysis  may  be  considerably  sim¬ 
plified  when  n  is  near  one.  For  the  fjjj  profile  we  write  Eqs.  (20)  and  (21)  in  the 
form 
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where  e  =  1  -  — .  Since  - = —  <  1  we  expand  the  binominal  term  in  the  powers 
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of  €  and  upon  integration  obtain,  to  first  order  in  e, 
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For  small  b,  Ij  and  Ig  may  be  written  as 
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As  indicated  previously,  the  Jj  and  J ^  integrals  are  duplicates  of  the  1^  and  I ^ 
integrals  with  b^  replaced  by  b.  We  then  have,  for  the  linear  f^  profile. 
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A  numerical  check  of  the  approximate  solutions  developed  in  this  section  were 
within  a  few  percent  of  their  more  exact  counterparts  in  Section  4,  even  for  e  as 
high  as  0.  75.  The  good  accuracy  obtained,  even  for  moderate  values  of  e,  is  due 
to  the  small  value  of  the  integrated  term.  It  is  estimated  that  the  first  neglected 
term  in  the  expansion  is,  in  general,  only  a  few  percent  of  the  first  (e  free)  term. 


6.  DETERMINATION  OF  REAL  HEIGHT 

The  calculated  values  of  h’(f)  are  presented  as  a  function  of  frequency  for  the 

2 

fj^  profiles  in  Figures  2  through  4  and  for  the  linear  f^  profiles  in  Figures  5 
through  7.  The  frequency  f  has  been  normalized  to  x  =  f/f  ,  where  f  is  the  max¬ 
imum  value  of  the  probing  frequency  f  for  a  given  curve.  The  parameter 
by  fv/fp  has  been  used  to  characterize  the  set  of  curves  for  a  given  figure.  The 
virtual  height  has  been  normalized  to  its  maximum  value  h'(fy).  Values  of  h'(f), 
aside  from  the  coefficients  a  and  g,  may  be  obtained  directly  from  Figures  8  and  9, 
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Figure  2.  Virtual  Height  h'(f)  for  the  fyg  Profiles  as  a  Function  of  Frequency 
x  for  Several  Values  of  Maximum  Frequency  Parameter  b;  Parameter  n  =  4. 
The  virtual  height  h'(f)  has  been  normalized  to  its  maximum  value  h'(f  ) 
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Figure  3.  Virtual  Height  h'(f)  for  the  fjyj  Profiles  as  a  Function  of  Frequency 
x  for  Several  Values  of  Maximum  Frequency  Parameter  b;  Parameter  n  =  2. 
The  virtual  height  h'(f)  has  been  normalized  to  its  maximum  value  h'(f  ) 


Figure  4.  Virtual  Height  h'(f)  for  the  f^j  Profiles  as  a  Function  of  Frequency 
x  for  Several  Values  of  Maximum  Frequency  Parameter  b;  Parameter  n  =  4/5 
The  virtual  height  h'(f)  has  been  normalized  to  its  maximum  value  h'(fv) 


Figure  5.  Virtual  Height  h'(f)  for  Linear  fw  Profiles  as  a  Function  of  Fre¬ 
quency  x  for  Several  Values  of  Maximum  Frequency  Parameter  b;  Parameter 
n  =  4.  The  virtual  height  h'(f)  has  been  normalized  to  its  maximum  value 
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Figure  6.  Virtual  Height  h'(f)  for  Linear  fN  Profiles  as  a  Function  of  Frequency 
x  for  Several  Values  of  Maximum  Frequency  Parameter  b;  Parameter  n  =  2. 

The  virtual  height  h'(f)  has  been  normalized  to  its  maximum  value  h'(f  ) 


Figure  7.  Virtual  Height  h’(f)  for  Linear  f„  Profiles  as  a  Function  of  Frequency 
x  for  Several  Values  of  Maximum  Frequency  Parameter  b;  Parameter  n  =  4/3. 
The  virtual  height  h'(f)  has  been  normalized  to  its  maximum  value  h’<f  ) 


which  give,  for  specific  values  of  n  and  b^,,  the  factor  h '(f  ),  by  which  the  normal¬ 
ized  values  of  h'(f)  are  multiplied.  The  a  and  g  coefficients  may  be  incorporated 
later  for  a  specific  probing  frequency. 

The  ionograms  chosen  for  determination  of  real  height  are  the  Boulder  and 


Garchy  ionograms  previously  utilized  in  Reference  3.  In  the  present  report  we 
omit  the  determination  of  initial  height  since  this  method  has  been  presented  in 
detail  in  Reference  3.  The  parameters  required  for  the  calculation  of  real  height 
are  obtained  from  a  comparison  of  a  given  ionogram  with  the  h'(f)  curves  developed 
herein.  The  h'(f)  curves  selected  for  the  Boulder  ionograms  are  plotted  in 
Figures  10  and  11,  and  those  for  the  Garchy  ionograms  in  Figures  12  and  13.  For 
greater  clarity  the  ionogram  data  are  plotted  as  discrete  points. 

As  an  example,  the  method  of  determining  the  real  height  h  is  presented  for 

the  Boulder  ionogram  at  the  higher  frequency  f  =  7.  65  MHz.  (For  convenience, 

values  of  h  and  h'  are  here  taken  relative  to  h  ): 

o 

From  Reference  3  h'  =  180  km 

h  =  76  km 

From  Figure  10  n  =  4,  linear  f^  profile,  b^  =  0.  90 

From  Figure  9  h'(fv)  =  0-557  . 

We  now  determine  the  parameter  a  for  the  parabolic  profile  by  comparing  the 
virtual  heights  for  the  ionogram  and  the  parabolic  profile: 

From  Reference  3,  Figure  3  a  -  0.  7 

From  Fqs.  (17)  and  (18)  g:  -  1.  667,  g2  =  1.  148 

h'  =  glg2ah'(fv) 

180  =  (1.  667)(1.  148)(0.  557)a 
a  =  168.  8  km  . 

For  the  parabolic  profile  the  real  height  is,  from  Eq.  (7), 

h  =  a[l  -  (1  -  bv)1/4] 
h  =  168.  8(0.  438) 

-73.9  km 


17 


0.2 


0.4 


0.6 


0.8 


1.0 


0 


X 

Figure  10.  Comparison  of  Boulder  Ionogram  for  f  =  7.  65  MHz  With  h'(f)  Curves 
Linear  Profiles,  n  =  4 


Figure  11.  Comparison  of  Boulder  Ionogram  for  f  =5.0  MHz  With  h'(f)  Curves; 
Linear  f^  Profiles,  n  =  2  ' 
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Figure  12.  Comparison  of  Garchy  Ionogram  for  f  =  4.  4  MHz  With  h'(f)  Curves; 
Linear  fN  Profiles,  n  =  4 


Figure  13.  Comparison  of  Garchy  Ionogram  for  f  =3.0  MHz  With  h'(f)  Curves 
Linear  f.T  Profiles,  n  =  2  v 


In  a  similar  manner,  the  following  results  were  obtained  for  the  Boulder  and 
Garchy  ionograms: 

Boulder,  f  =  7.  65  MHz,  a  =  0.  7,  by  =  0.  9,  n  =  4  (linear  f^) 

h(calc. )  =  73.9  km,  h(Ref.  3)  =  76  km 

Boulder,  f  =  5.  0  MHz,  a  =  0.  6,  bv  =  0.  75,  n  =  2  (linear  f^) 

h(calc. )  =  34.  1  km,  h(Ref.  3)  =  33  km 

Garchy  f  =  4.  4  MHz,  a  =  0.  6,  by  =  0.  99,  n  =  4  (linear  f^) 

h(calc.)  =  92.0  km,  h(Ref.  3)  =  98  km 

Garchy,  f  =3.0  MHz,  o  =  0.  5,  b  =  0.95,  n  =  2  (linear  fM) 

J  v  v  N 

h(calc. )  =  27.9  km,  h(Ref.3)  =  31  km 

As  seen  from  Figures  10  through  13,  the  ionograms  could  be  fitted  reasonably 
well  to  the  parabolic  profiles,  except  for  the  Garchy  ionogram  at  the  lower  fre¬ 
quency  of  3.  0  MHz.  For  this  case,  the  ionogram  bends  too  rapidly  at  the  lower 
values  of  h'  before  following  the  trend  of  the  parabolic  profiles.  However,  other 
trial  fits  of  the  ionograms  to  the  parabolic  profiles,  yielded  real  heights  fairly 
close  to  the  values  obtained  herein,  indicating  that  a  reasonable  fit  should  yield  a 
representative  value  of  the  true  height.  Improvements  in  the  fitting  of  the  iono¬ 
grams  to  the  parabolic  profiles  and  in  the  accuracy  of  the  real  heights  can  pre¬ 
sumably  be  achieved  by  calculating  for  several  intermediate  powers  of  f^  and  for  a 
greater  range  of  n. 

7.  CONCLUDING  REMARKS 

A  direct  method  for  obtaining  real  heights  from  virtual  height  data  has  been 
obtained  by  developing  a  series  of  generalized  parabolic  profiles  with  which  a  given 
ionogram  may  be  compared.  Ionograms  from  Boulder,  Colorado,  and  Garchy, 
France,  were  chosen  for  the  calculations.  The  results  obtained  were  in  fairly 
good  agreement  with  those  obtained  from  previously  reported  accurate  numerical 
calculations. 

The  present  method  affords  a  generalization  of  the  procedure  presented  in 
Reference  3,  in  which  virtual  height  data  was  represented  by  a  parameterized 
series  of  polynomials  from  which  real  heights  could  be  obtained  by  the  evaluation 
of  simple  integrals.  Since  the  parabolic  profile  has  a  vertical  slope  at  its  maximum 
electron  density,  it  affords  a  more  realistic  fit  to  ionospheric  electron  density  data 


